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ABSTRACT 



The omni-present low frequency wave motion (30-300 sec) contains a substantial 
fraction of the total wave energy inside the surf zone. A more complete description of 
nearshore wave processes considers incident short period wind waves superposed on, 
and interacting with, long standing waves. The wind waves are modulated in 
amplitude, wavenumber and direction due to relatively slowly varying depth changes 
caused by the long waves. The energy in the wind wave band is enhanced by side band 
growth at the sum and difference frequencies of short and long waves (order 15% at 
the shoreline). The modulation is identified in the analysis of field data as a positive 
correlation between the long waves and the wind wave envelope near the shoreline. 
Considering oblique incident waves, a steady longshore current showing a 
non-vanishing current at the shoreline is found as a result of the non-linear interaction 
between monochromatic incident and infragravity waves. An analytical solution 
describing the unsteadiness of the longshore current is developed. Applying the derived 
longshore current solution, longshore sediment transport is reformulated to include the 
infragravity waves, giving improved comparisons with field measurements. 
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I. INTRODUCTION 



Sea waves are generated by wind, and as they propagate towards the shore they 
deform due to shoaling and may refract if they arc oblique to the bottom contours or if 
they meet a current. The waves grow in height until they become unstable and break, 
dissipating their energy. The processes of breaking waves on a sloping bottom make 
the surf zone an extremely dynamic area where different observable phcnomenac take 
place, such as longshore current, sediment transport and wave runup (setup and 
swash). Due to this complexity, the proper modeling of wave motion in the nearshore 
zone has been the goal of many investigators, and yet much remains to be done before 
a satisfactory understanding is achieved. 

The surf zone can be characterized dynamically by three regions: the outcr.thc 
inner, and the runup region (Svendson et al., 1978), (see Fig 1.1). The outer region 
follows the breaking point and is dominated by incident breaking waves. The breaking 
waves are of either the plunging or spilling type, and are characterized by a rapid 
transition of wave shape with a horizontal surface roller. In the inner region, the wave 
forms break down into smaller scales of random and turbulent nature and eventually 
become similar to moving bores, or hydraulic jumps. The run-up region is 
characterized by wave setup and swash oscillation running up and down the beach 
face. In all regions, the wave height and the water motion are strongly locally 
controlled by the depth. The depth controlled waves are referred to as "saturated." 

In general, surf zone dynamics are found to differ dramatically between reflective 
(steep) beaches and dissipative (gentle slope) beaches (Bradshaw, 1980). Typical 
features of the reflective beaches are high reflection of incident waves, collapsing or 
surging breakers on the beach face, and subharmonic resonances (c.g. Guza and Davis, 
1974; Guza and Inman, 1975; Huntley and Bowen, 1975b, 1978; Wright et al., 1979). 
Dissipative beaches arc dominated by spilling breakers, pronounced low frequency 
oscillations and the presence of multiple parallel longshore bars (c.g. Huntley and 
Bowen, 1975a; Short, 1975; Sasaki and Horikawa, 1975; Holman et al., 1 9 7 S ; Symonds 
et al., 1982). This study will focus on the dynamics of dissipative beaches in which the 
inner region is assumed to be a transition region between regions dominated by 
incident breaking waves and low frequency swash. 




Figure 1.1 Surf zone regions and coordinate system. 
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In the past decade, considerable attention in the literature of beach and 
nearshore dynamics has been devoted to the kinematics of low frequency waves. These 
waves, which have periods ranging from 30-300 sec, are referred to as "infragravity 
waves." The observations show that the dominant wave motion closer to the shore is 
not normally at the incident wave frequency band. It is also found that the low 
frequency wave energy in the very nearshore region contains a substantial fraction of 
the total wave energy (Inman, 196S a,b; Suhayda, 1974; Goda, 1975; Huntley and 
Bowen, 1975b; Sasaki and Ilorikawa, 1975,1978; Huntley, 1976; Wright et ah, 
1978, 1979, 19S2; Bradshaw 19S0; Holman, 1981; Huntley et ah, 1981; Guza and 
Thornton, 1982,1 985a, 19S5b; and others), (see Fig 1.2 for example). The valley in the 
energy spectra, where the energy level reaches minimum, suggests an empirical cutolf 
frequency around 0.03 Hz that separates the two wave bands. Therefore, a complete 
description of the wave field in the surf zone should include two relatively separated 
bands of wave frequencies : (1) incident sea-swell waves (short period waves) with 
periods of 1-30 sec, and (2) infragravity waves (long period waves) with periods of 
30-300 sec, which no longer can be neglected inside the surf zone. 

Prior to the 1950's, the incident wave characteristics were the main concern when 
measuring and analyzing waves in the nearshore zone. Munk (1949) and Tucker (1950) 
were the first to point out the existence and significance of low frequency motion and 
called the phenomenon "surf beat." More recently, with renewed interest in nearshore 
processes and with major improvements in the instrumentation, it has become apparent 
that the low frequency wave energy can be as important as the incident wave energy 
when approaching the shoreline. Another motivation for the study of the infragravily 
waves is the coupling between the infragravity length scales (order of hundreds of 
meters) and the length scales associated with most of the rhythmic beach features such 
as sand bars, beach cusps and crescentic bars. Many investigators have related the 
beach morphology changes to the existence of low frequency motions by which a 
mechanism of forming and maintaining such topographic features is provided (Bowen 
and Inman, 1971; Guza and Inman, 1975; Bowen, 1980; Holman and Bowen, 1982: 
Bowen and Huntley, 19S4; and others). 

In analyzing surf zone data, there are often anomalous patterns of water motion 
closer to the shoreline when compared with present short wave models. For example, 
the spatial distribution of wave height measured in both the laboratory and field often 
shows unexpected wave height at the shoreline (e.g. Suhayda and Pettigrew, 1977: 
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(Das/ 2 uo) A 1 ISM 30 A9U3N3 



Figure 1.2 Cross-shore velocity spectra on 20 Nov.. 1978, Torrey Pines 
(after Gu/.u and Thornton, 1985b). 
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FREQUENCY (HZ) 



Battjes and Stive, 1985), (see for example Fig 1.3). Bowen, ct al. (1968), in measuring 
the setup, found a residual wave height always present at the shoreline and modified 
their setup formula at the beach face to fit the data. Suhayda and Pettigrew (1977) 
measured the average wave velocity across the nearshore zone and found it to decrease 
shoreward until the sea-swash limit of the run-up zone was reached, and then to 
increase onshore, supporting Waddell's measurement (1973). Thus, the wave height 
and celerity do not vanish at the vicinity of the shoreline, invalidating the widely 
applied short wave saturation assumption. 

The assumptions of wave saturation and steady state conditions for both 
monochromatic and random waves, results in a zero current at the shoreline. In 
addition, the steady longshore current assumption is often unrealistic, especially with 
the presence of irregular oscillating waves (Wood and Meadows, 1975). This is often 
seen in field data that show longshore velocity oscillations at longer wave periods than 
for incident waves (c.g. Holman and Bowen, 1984; Guza and Thornton, 1985a; 
Oltman-Shay and Guza, 1986). Using an empirical orthogonal function (FOF) analysis 
on 64 data sets covering a wide variety of wave conditions, Guza et al. (1986) obtained 
objectively best fit to the longshore current distribution. The analysis showed the 
"classical parabolic" shape for the longshore current distribution (Bowen, 1969) over 
most of the surf zone, but with significant residual velocity at the shoreline (Fig. 1.4). 
In the analysis of the same data, Wu et al., (1985) found a poor agreement between 
longshore current theory with velocity going to zero at the shoreline and data closer 
to the shoreline. All present models predict zero current at the shoreline, and no model 
is available to explain the unsteadiness of longshore current. 

Improving the longshore current description is important because it provides the 
driving means of transporting sediment alongshore (Longuet-lliggins, 1972). The 
available longshore sediment transport models (c.g. Bijkcr, 1971; Thornton, 1973; 
Komar, 1977; Sawaragi and Dcguchi, 1979; Bailard, 1981), all show a longshore 
transport maximum closer to the breakcrlinc, which agrees with the longshore current 
distribution. The shortcoming of these models is their failure to predict any moving 
sand at the shoreline regardless of the wave conditions. Zcnkovich (1960) measured 
the suspended sediment transport distribution across the surf zone using fluorescent 
sand tracers. Sand transport was found greater over bars where energy dissipation is a 
maximum due to wave breaking and also at the shoreline, (sec Fig 1.5). Measuring 
longshore sand transport distribution using sand tracers, Kraus ct al. (1981) found a 
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I'ieurc 1.3 Root mean square wave height H vs offshore distance x. 
Solid curve represents a short wave modcTand tmrstars arc laboratory data 
showing non-zero wave height at the shoreline (after Uattjes and Stive, 1 985). 
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distinct bimodal distribution having maxima at the breakerline and near the shoreline. 
The same results were obtained by Downing (1984) when measuring the suspended 
longshore sediment distribution. Abdelrahman (1983) compared longshore sediment 
transport at Santa Barbara beach with his predicted longshore sediment transport 
model based on the widely accepted energetic concept by Bagnold (1963). The model 
did not predict the large transport in the vicinity of the shoreline since theoretical short 
wave description gives zero amplitude at the shoreline. It was suggested a primary 
cause of these differences was the presence of infragravity waves, which have not been 
taken into consideration before. 

In recognition of the importance of the infragravity waves in the surf zone, an 
obvious next improvement in the description of waves kinematics is to combine the 
infragravity waves with the incident sea-swell. Asa result, the incident waves would be 
modulated in amplitude, wavenumber and direction. The purpose of this dissertation is 
to improve the understanding of nearshore processes by generating and applying such a 
composite wave description. 

As will be seen, the interaction of the short waves with the long waves results in 
energy not only at their respective frequencies, but energy transferred to side bands of 
the incident waves at the sum and difference frequencies. It is the energy and 
momentum at the side bands that is primarily responsible for the inner surf zone 
dynamics. The dynamical changes are demonstrated by the successful treatment of 
using the first order wave description in the momentum flux equation to obtain a 
second order dynamical effect. However, inclusion of infragravity waves does not 
mean introducing more independent parameters since the infragravity wave energy is 
assumed related to the incident wave field. 

In the following chapter, a description for the short and long period waves is 
presented. The short period waves are described by linear (Airy) theory, while the long 
period waves on a sloping bottom are described by either cross-shore standing waves 
(2-D) or edge waves (3-D). In chapter 3, the changes in the amplitude and 
wavenumber of the short waves superposed on long waves are discussed. An analytical 
expression describing the modulated surface elevation of the short waves is developed 
along with their energy and momentum. Changes in the wave kinematics are 
considered with the long period waves treated as a slowly varying current. In chapter 
4, narrow banded wave data are analyzed to examine the hypothesis of short wave 
modulation by the long waves. The dynamical effects of considering the infragravity 
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waves inside the surf zone are investigated. Analytical models for calculating wave 
setup, longshore current and sediment transport are developed in chapter 5. An 
analytical solution to the unsteady longshore current is derived. Model validation and 
comparison with other models are included in chapter 5. The conclusions of this study 
are presented in chapter 6. 
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Figure 1.5 Distribution of longshore velocity and sediment transport 
across the surf zone (al ter Zenkovich, 1960). 
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II. WAVES IN THE NEARSHORE ZONE 



A wave record of incoming waves usually contains short period waves (1-30 see) 
and long waves (30-300 sec). To study first order modulations of short waves 
superimposed on long waves, the principle of superposition is applied. Therefore, the 
non-linear wave theories are excluded in this study. A chosen Cartesian coordinate- 
system is illustrated in Fig. 1.1 where the positive X-axis is directed offshore, the 
positive Z-axis is vertically upward and Y is alongshore. This coordinate system is 
convenient for describing the offshore decay of the standing long waves. Linear 
descriptions of both sea-swell and infragravity waves are discussed separately. Short 
waves are derived from linear wave theory over a flat bottom, then using shallow water 
approximation to describe short waves in the nearshore zone. Infragravity waves are 
derived from shallow water equations on a sloping bottom. 



A. SHORT PERIOD WAVES 

In general, waves in a viscous fluid propagate over irregular topography of 
varying permeability. In most eases, the main body of the fluid motion is reasonably 
assumed irrotational since the viscous effects are usually limited to a thin boundary 
layer near the surface and the bottom. Surface waves arc considered a boundary value 
problem. 

Linear wave theory (small amplitude waves) is considered a first order 
approximation to the theoretical description of wave behaviour. It is assumed that the 
wave amplitude , a, is very small compared to both the wavelength L and the local 
water depth h, (a < < L and a < < h). 'flic fluid is assumed incompressible, 
homogenous, in viscid and also irrotational, by which a velocity potential (p should exist 
to represent the field of the flow: 



u 



7 $ 



( 2 . 1 ) 



where u is the velocity vector, liquation 2.1 satisfies the continuity equation V ,u = 0, 
by which the governing differential equation (Laplace equation) is obtained: 



V 



2 



c 



| ~h <_ z £ n 
1 -<« < x < 00 



( 2 . 2 ) 



With the proper kincinatical and dynamical free surface and rigid flat bottom boundary 
conditions, the method of separation of variables is used to solve equation 2.2 Then, 
the velocity potential is given by 



4 > 



acj cosh — k ( h + zj _ cos (^ x + u t) 
m sinh kh 



(2.3) 



where a is the wave amplitude, 0) = (2it/T) is the radial frcqucncy, and k = (2rc/L) is 
the wavenumber. 'I he surface elevation t| is periodic in both space and time. 



n = a cos (k^x^ + ut) , i = 1,2 



(2.4) 



where i = 1,2 refers to conditions in the X and Y-dircctions. The description of the 
wave motion is completed by the dispersion relationship, 



2 

CO 



gk tanh kh 



(2.5) 



which arises from the dynamical free surface boundary conditions. The horizontal 
velocity under the wave is given by 



u 



cosh k (h+z ) 
aj sinh kh 



cos 



(kx + cat) 



( 2 . 6 ) 
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and the vertical velocity 



W 



aw 



sinh k (h + z ) 
sinh kh 



si n ( kx + wt ) 



(2.7) 



Swell waves propagating from deep water are nearly sinusoidal with long low 
crests. As waves approach shallow water, they start to feel the bottom and deform at 
relative depth h/L < 1/2. The wavelength decreases while the wave height generally 
increases, and the wave period remains invariant. The hyperbolic functions have 
convenient deep and shallow water asymptotes by which an approximate solution can 
be obtained. In shallow water, the linear gravity waves arc non-dispersivc since the 
dispersion relationship reduces to 



2 

Cd 



g k 2 h 



( 2 . 8 ) 



and the horizontal and vertical velocities simplify to 



u = a /g/h cos (kx + ut) 



(2.9) 



W = 



ao ( 1 




sin (kx + ut) 



(2.10) 



On dissipative beaches, long crested swell waves break when the wave steepness 
increases and the velocity of the water particles at the crest exceeds the phase speed of 
the wave form. In the inner surf zone, bores retain their relatively long crested form 
and progressively decrease in height as they advance onshore. Due to turbulence and 
non-linearities associated with the breaking process, the motion of the water particles 
after breaking can no longer be described analytically. As a first approximation, the 
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wave or bore height, II, at any location inside the surf zone is limited by the local 
depth according to the saturation assumption, i.c. 



H = 2yh 



( 2 . 11 ) 



where y is a constant of order unity. This relation holds for spilling breakers, the most 
common breaker type on natural beaches. The breaker index y is important in 
describing the breaking waves used in longshore current and sediment transport 
models. This linear relationship is expected on dimensional grounds (Longuct-1 liggins, 
1972) and gives a reasonable description of laboratory results for monochromatic 
waves (c.g. Galvin and I;agleson,1965) and field measurements during conditions of 
both broad and narrow spectral distributions of wave energy (Thornton and Guza, 
1982), (Fig. 2.1). 

Very close to the shoreline in the run-up region, the fluctuating component of the 
wave runup about the mean water level is known as "swash oscillation." Michc (1951) 
hypothesized that monochromatic waves within the surf zone may be composed of 
both a progressive component, which is dissipated onshore and has zero shoreline 
amplitude, and a standing component, which has its maximum amplitude at the 
shoreline. Many laboratory studies (Moraes, 1970; Battjcs, 1974; Guza and Bowen, 
1976; Van Dorn, 1976) report that the shoreline swash amplitude a $ is limited by the 
swash parameter c s , 



£ 



S 




( 2 . 12 ) 



where tan (1 is the beach slope and c $ varies from one experiment to the other in the 
range from 1 to 3. Guza and Thornton (1982) attributed this variation to the 
ill-delincd nature of the backwash (run-down). Hunt (1959) related the total runup R v 
, defined as the sum of setup and half the swash amplitude, to the surf similarity £.q 
with an empirical constant M: 
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Figure 2.1 Wave heights inside the surf zone defining breaker index 
and showing saturation (after Thornton and Gu/.a, 1982). 
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(2.13) 



R V /H sig - M e 0 - M tan p / V (H 0 / L 0) 



where H.:_ is the significant wave height and the subscript 0 denotes the deep water 
conditions. 

Huntley et al. (1977) showed a swash "saturation" condition at incident wave 
frequency, i.e. with increasing wave height, the steady set-up will increase but the 
swash amplitude will not. Guza and Thornton (1982) found a spectral decay at 
incident wave frequencies, indicating that energy levels are independent of incident 
wave height, although Huntley et al. (1977) have shown co*^ dependency. The swash 
saturation condition is found to be true for both monochromatic and random swash 
fluctuations (Guza et ah, 1984). Swash at low frequency is found to be unsaturated, 
i.e. the swash amplitude increases with increasing incident wave height. Most of the 
low frequency energy is found to have infragravity wave periods. Thus swash spectra 
show a saturated region at incident wave frequency and an unsaturated region at low 
frequencies that dominates the whole spectrum. 

To evaluate the amount of runup, swash studies have been mainly concerned 
with measuring the swash amplitude at the shoreline without putting much emphasis 
on the generation and types of the dominant low frequency motions. It is of interest 
now to investigate the forms of this low frequency water motion that give maximum 
amplitude at the shoreline. 



B. LONG PERIOD WAVES 

The long period waves presented in this work are limited to waves having a time 
scale of 30-300 seconds where gravity is the primary restoring force. 

1. Infragravity Wave Models 

The principal wave types that may contribute to the low frequency motion on 

beaches can be summarized by their mathematical formulation as follows: 

(i) 2-D model : the wave motion is mainly in the cross-shore direction and is 

known as surf beat or "leaky modes." 

(i) 3-D model : the wave motion has longshore variation and is known as edge 

waves or "trapped modes". 
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The mathematical forms of both waves arc derived from the linear shallow 
water wave equation on a sloping bottom 




- 9% (h 



3 0 ^ 
3x J 





0 



(2.14) 



In the 2-D model , the wave energy is reflected offshore and the solution 
obtained is a standing wave expressed in terms of the zero order Bessel function of the 
first kind, .1 q , (Lamb, 1932; Frcdriehs, 194S) 



$(x,t) = Jq (X ) sin wt 



(2.15) 



9 1/9 

where y = (4co x/ g tan p) ' *■ and a is the amplitude at the shoreline. It is noted that 
the Bessel function of the second kind, Yq , is excluded since it is unbounded as waves 
approach the shoreline. The surface elevation and the associated orbital velocities arc 
given by 



n(x,t) = a j (>:)cos wt 



(2.16) 



u(x,t) = a /g/h J^(X) sinwt 



(2.17) 



W (x , t ) = - au Jq (X) sin ut 



(2.18) 



In the 3-D models, the edge waves are free waves that propagate along the 
coast with their energy trapped onshore by refraction, i.e. they do not radiate energy 



offshore. Stokes (1846) was the first to provide an analytical solution to describe edge 
waves using a small amplitude wave assumption. Eckart (1951) solved equation 2.14 
analytically in terms of Eagucrre polynomials E n (2ky x). 'flic solution of progressive 
edge waves is 



= ®SL 



-k x 

, y 






L n ( 2 ky x) cos (k^y 



- uit ) 
n 



(2.19) 



where ky is the longshore wavenumber, and n is an integer number which represents 
the edge wave mode number. Each discrete mode contains finite energy, and the n-th 
mode edge wave will have n zero-crossings after which the amplitude decays 
exponentially. Edge waves obtained by Eckart (1951) must satisfy the dispersion 
relationship, 
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CJ 



n 



q k ( 2n+l ) tan 3 

^ y 



(2.20) 



which requires (2n + I)tanp < < 1 to insure that the solution achieves its limiting value 
as ky x oo while still in shallow water. 

Ursell (1952) was able to obtain a set of "exact" edge wave solutions on a 
sloping bottom using the small amplitude wave theory without recourse to the shallow 
water approximation. His solution satisfies a slightly different dispersion relationship 
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n 



q k sin ( 2n+l ) 6 

J y 



( 2 . 21 ) 



and the solution lies within the range ((2n+ 1)|) ) ^ 71/2 to insure a cutoff mode for a 
given beach slope. Ursell (1952) showed that Stokes'( 1 846) edge wave solution is only 
the zero mode of his solution. For gentle beach slopes and low mode numbers, both 
dispersion relationships 2.20 and 2.21 arc in agreement, suggesting that the shallow 
water approximation is appropriate under these conditions, (Guza and Davis, 1974). It 
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is noted that the approximate Eckart (1951) solution is more commonly used in 
comparisons with field data due to its simpler mathematical formulation. In the field 
data analysis, a difficulty in discriminating between leaky and edge wave modes for 
small offshore distances arises because of the similarity of their cross-shore structures. 
Also, both waves have maximum amplitude at the shoreline. Fig. 2.2 shows such 
similarities between higher mode edge waves and reflected waves. 

In summary, the low frequency motions may be classified as edge waves, leaky 
modes, forced waves or a combination as shown in Fig 2.3, (after Symonds, 19S2). A 
complete set of free wave solutions of discrete edge waves occurs for CO" < gk y while a 
continuum of leaky mode exists for > gky . It is noted that the forced waves may 
lie between the edge wave modes, but they do not satisfy the dispersion relationship 
and arc weaker than the free modes since edge waves may grow rcsononantly as they 
are excited. 

2. Generation of Infragravity waves 

Munk (1949) and Tucker (1950) suggested that the long period waves may be 
caused by an excess of mass carried forward by groups of high swell. The swell waves 
are assumed to be destroyed on the beach, and the excess mass transport is reflected 
back as a free long wave. They tested this hypothesis by comparing offshore pressure 
records, delayed by the calculated travel time of the incident wave envelope. Tucker 
(1950) observed a negative maximum correlation between incident wind wave groups 
and the infragravity waves at a time lag equivalent to the travel time of a forced wave 
moving with the wind wave group velocity and traveling back upon reflection as a free 
wave. 

Longuct-I liggins and Stewart (1962) showed that the variations in the 
radiation stress due to incident wave grouping drive a forced wave such that a 
depression of the mean water level (MWL) will occur under incident high waves, and a 
corresponding rise in MWL occurs under low waves. This forced oscillation of second 
order is associated with groups until the incident waves break onshore and release a 
free wave offshore, which theoretically supports the above observation. 

Gallagher (1971) suggested edge wave excitation through a non-linear 
interaction whereby two incident waves can transfer energy into an edge wave if the 
difference frequency and longshore wavenumber of the two interacting wind waves 
satisfy the edge wave dispersion relationship, equation 2.21 . This resonance has been 
observed in the laboratory (Bowen and Guza, 1978). 
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Figure 2.2 Modes of edge waves and reflected waves, <p is the velocity potential. 
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Incident waves may transfer energy to either leaky or trapped modes through 
many possible interactions. One of these is given by Guza and Davis (1974) where 
energy is transferred into edge wave modes through a non-linear interaction involving a 
normally incident standing wave on a sloping bottom. Bowen and Guza (1978) showed 
energy transfered into leaky modes (to - > gk v ) by edge waves. 

Another mechanism for low frequency wave generation is suggested by 
collision between the stronger component of backwash and the shoreward moving 
bores that results in roll waves which last for few seconds. On analyzing swash data, 
Waddell (1973) found that the collision between uprush and backrush plays an 
important role in reducing the extent of swash runup waves. Furthermore, he 
suggested that water percolation into the beach serves as a low-pass filter alfccting the 
swash frequency. Bradshaw (1980) suggested bore-bore capture as a mechanism for 
long wave generation in the swash zone. Mase and hvagaki (1984) showed a 
considerable shoreward decrease in the ratio between the number of runup waves to 
the incident waves. 

Symonds et al. (1982) developed a 2-D model for generating long waves which 
is forced by a time-varying break-point caused by the groupness of the incident wave 
field. The generated free long wave at the group frequency has a standing wave 
structure shoreward of the breakpoint and an outgoing progressive wave structure 
seaward of the breakpoint. 

In the above survey, both short and long wave models are described separately 
showing different characteristics. It is of interest now to combine the two waves and 
study the changes in the incident wave (short wave) field due to the presence of longer 
period waves. 
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Figure 2.3 Dispersion relations for infragravity waves. 
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III. SNORT WAVE MODULATION DUE TO LONG WAVES 



When short period waves ride on much longer waves, they arc modulated in 
amplitude, wavelength and direction. Their amplitude tends to increase at the crests of 
the long waves and decrease at the troughs. The changes in wavenumber and direction 
depend mainly on kinematical considerations. In this chapter, the changes in 
wavelength, direction and amplitude arc investigated. A development lor the 
corresponding changes in energy and momentum flux then follows. 

A. CHANGES IN THE WAVE AMPLITUDE 

The amplitude modulation of a short progressive wave riding on a much longer 
wave was pointed out by Unna (1942, 1947). Unna (1947) showed a contraction of the 
wavelength and an increase in the amplitude of short waves at the long wave crests. 

Longuct-1 liggins and Stewart (1960) used a perturbation analysis to examine the 
non-linear interaction between short and long progressive waves over a horizontal 
bottom. They described the surface elevation 



n = a s sin i|/ g + a £ sin ip z 



(3.1) 



where »| / = (kx- to t + 0 ), is the phase function and 0 the phase shift. The subscripts 
s and £ refer to the short and long wave respectively. Using Stokes' method of 
approximation to the second order, they derived a general form for the modulated 
surface elevation 



' = a (1 +P)sin + a Qcos 

S b ^3 



(3.2) 



where P = a^ k^ sin \\t^ and Q = -a^k s cos . Using the first term in the right hand 
side of equation 3.2, Longuet-I liggins and Stewart expressed the modulated short wave 
amplitude as 
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(3.3) 



= a 



{1 + 



a i k s 



j coth k ( h 



x tanh k „h) sin y. ) 

4 *. ** 



and a modulation in the short wave wavenumber, k' , as 



k’ - k { 1 + a^k^ cothk^h simj^j 



(3.4) 



In shallow water, equations 3.3 and 3.4 reduce to 

a ’ = a s 11 + ! -k sin,|, i ) 



k’ = k g { 1 + -^sintf £ } (3.6) 



These amplitude and wavenumber modulations arc explained by the work done by the 
long waves against the momentum of the short waves. This work is converted into 
short wave energy and therefore produces a steepening of the short wave at the crests 
of the long waves, i.e., the energy is redistributed along the wavelength of the long 
wave. Longuct-I Iiggins and Stewart (I960) came up with equations similar to 3.3 and 
3.4 for a progressive short wave superposed on a long period standing wave in which ai 
was described as 2aj . 

I. Determination of the Modulated Wave Profile 

It is of interest now to show that the refractive changes in amplitude and 
wavenumber due to the non-linear interaction may result in a modulated short wave 
expressed as an infinite sum of components at different frequencies. Then, the 
properties of the modulated wave are investigated to provide a better understanding of 
the modulation process. The procedure is to let the modulated short wave, designated 
bv primes, propagate in the same wave direction and satisfy 

n’ = a’ sin (k'x - w ' t ) (3.7) 
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Substituting equations 3.3 and 3.4 into equation 3.7, and letting u) = w s , yields 



n ' = a' sin (0 + sin 4^) 



(3.8) 



where = a^ k s x coth k£ h is defined as a non-dimensional wavenumber 
(wavenumber modulation factor). liquation 3.8 is a general form of amplitude and 
wavenumber modulation. liquation 3.S reduces to the simple form rj = a sin vj/ 
when nij. < < 1. Expanding equation 3.8 using trigonometric identities, gives 



n ' = a'{sin cos (m^ sin + cos sin (m^ sin tj^) } (3.9) 



Using Fourier series expansions to express the sinusoidal arguments (Abramowitz and 
Stegum, 1964), gives 



n' = a' l J n < m k.) sin ( + nij^) (3.10) 

n=-°° 

In the above analysis no assumptions have been made regarding either applicable 
regions or the initial properties of the two interacting waves. Equation 3.10 offers a 
means of expressing the modulated short waves as a combination of spectral 
components based on monochromatic input. The short wave is viewed now as a 
carrier subjected to a simultaneous amplitude and phase modulation, (Panter, 1965). 
Changing the limits of summation, equation 3.10 becomes 



n' = a'J 0 (m k )sin ij> g 



(3.11) 
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The first term has the form of the solution for waves at the incident wave frequency 
propagating on a sloping bottom but with a different Bessel argument (sec equation 
2.16 ). The remaining terms represent the side bands at sum and difference frequencies 
(\|/ s ± n\| ij^) where they arc phase locked. This solution is dependent on the argument 
m^ and converges rapidly as n increases. 

In shallow water, the Bessel function argument nij, reduces to 2 (a^ /h) k $ x 

for the ease of progressive short waves superposed on long standing waves, 'flic 

argument nij. is found to increase onshore (Fig 3.1), where the computed values of m^ 

arc evaluated until closer to the shore (x = 1 cm) before reaches oo . Due to the 

onshore increase of m^ , ,l n (m^ ) decreases onshore, oscillating rapidly and reaches 

zero at the shoreline (Fig 3.2). Note the difference in behaviour of the Bessel functions 

in equation 3.11 and 2.15, where the Bessel argument allows the amplitude of the long 

standing wave to reach a maximum at the shoreline. For breaking short waves on a 

long period standing wave, the depth limited breakers is utilized to describe a s (= y h) 

in equation 3.5 . Substituting a' into 3.1 1, yields 

3 

n’ = YhlJ 0 (r\)sin 4> s + j -%■ ^(n^lcos iftj 

+ Yh[J n (m k ){sin(;p s +n^) - (-l) n sin (if> s -rty £ )} 

(3.12) 

' S -F (J n-l V - J n+1 < V ){oos<l| 's +n ‘ l 'n) " 'V'V 11 

Two distinct features of this solution arc demonstrated in the amplitude squared 
spectra (Fig 3.3). First, symmetric side bands arc generated about the carrier (Fig 3.3 
a), 'flic second feature is the transfer of short wave energy to the upper and lower side 
bands as the depth decreases (Fig 3.3 b), which causes broadening of the wave spectra. 
It is also noted that the amplitude of the spectral components decrease shoreward. 
The instabilities associated with the side bands generation arc attributed to their spatial 
structure and phase variation. It is shown that the side bands energy is modified by 
the decreasing depth. 'I he modulation may contribute to the infragravity waves by 
strengthening the low frequency energy component closer to the shoreline, which has 
previously been observed. Waddell (1973) presents evidence of nonlinear transfer of 
energy from high to low frequency due to collision, where successive bores begin to 
overrun each other. However, Huntley and Bowen (1975b) suggested that energy 
transfers arc due to the interaction of many irregular spaced breakers in the surf zone. 
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Figure 3.1 Wavenumber modulation factor mi, vs non-dimensional distance x/x h , 

where x^ is the surf zone width. D 
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Sawaragi and Ivvata (1974) studied experimentally the wave transformation after 
breaking and observed energy transfer from the carrier into higher frequeney waves by 
an unknown mechanism. The present analysis may be used to illustrate some of these 
findings, liquation 3.11 shows that the beach docs not eliminate the incident wave 
frequency and that the modulated wave energy components and the incident wave 
frequency vary spatially. 



B. CHANGES IN THE WAVE KINEMATICS 
1. Background on Wave Kinematics 

At any instant of time, a wave front is defined by the phase function vj/ = 
constant, which is the equation of a family of parallel planes with normal vector k. As 
time increases, these planes move with the phase speed in the direction k. 

In general, the dispersion of the wave motion may be written 



cj = to ( k , h ) 



(3.13) 



where the frequency varies with both the wavenumber and the local water depth. The 
local wavenumber and frequency arc defined by the gradient of the slowly varying 
phase function in space and time: 



& = Vip 



(3.14) 



oj = 




(3.15) 



Since the curl of a gradient is zero, it follows immediately from equation 3.14 that the 
local wavenumber vector in space is irrotational, i.c., 
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Figure 3.2 Behaviour of the Bessel functions J () (nu\ 
(solid curve) and Jj (mj.^ (dashed curve) across the suifkbne. 
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0 



(3.16) 



V * k 



Cross differentiating equations 3.14 and 3.15 to eliminate vj/ yields 



Equation 3.17 is known as the kincmatical conservation of wave density, where the rate 
of change of wavenumber is balanced by the convergence of the frequency. Sometimes 
3.17 is also called the "conservation of crests equation," since the crests arc neither 
created nor destroyed and their total number must be conserved. 

Rewriting equation 3.17 in tensor notation and expanding the convergence 
term using the dispersion relationship in equation 3.13 gives 




D i 



0 , i = 1 , 2 



(3.1S) 



Equation 3.18 can be rewritten as 




0 , i = 1,2 



(3.19) 



where d kj jd Xj = d kj jd Xj from equation 3.16 , and the group velocity is defined by 



C 




(3.20) 



40 




"o 




Figure 3.3 Amplitude square spectra of the modulated short wave 
at offshore distance: (ajx = 40.0 m (b)x = 0.5 m. 
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2. Kinematics of waves on current 

In studying the problem of short gravity waves superposed on waves of a 
much longer wavelength and period, the long wave may be modeled as a slowly varying 
current U(x,t), (e.g. Longuct-1 liggins and Stewart, 1960, 1961; Garrett and Smith, 
1976). A coordinate system that moves with the current velocity is chosen, by which 
the observed frequency to , passing a fixed point, is given by 



to ' = to + k • U 

s s x 



(3.21) 



where 0) $ is the intrinsic frequency, (Ursell, 1960; Whitham, 1960). The subscripts s 
and 1 refer to short and long wave properties. The second term on the right hand side 
of 3.21 is due to the relative motion of short waves with respect to the chosen 
coordinate and is known as the doppler shift effect. 

T he short wave will be modulated by the presence of the current, and 
therefore the modulated short wave (denoted by primes) is assumed to satisfy equation 
3.17, i.c., 

+ Vu’ = 0 (3.22) 

a 



where k' , the observed (modulated) wavenumber, is unknown and assumed unsteady. 
The observed frequency to is given by equation 3.21 . Expanding equation 3.22 gives 



3k! 3k 

+ C — 

3 1 g. 3x. 

3 



. 3k. U 

3u) 3 n ^ i 

3h 3x^ 3x^ 



0 



(3.23) 



where the quantities kj , eg and to arc for the short waves, and LJ is the current (long 
waves). 

Assuming a steady intrinsic short wave train; i.e., d k s / d t = 0, equation 3.23 
may be simplified by the use of equation 3.19 to give 



42 



0 



(3.24) 




o :< . 
l 



which describes the kinematic conservation of the modulated wave density. 

3. Applications 

'flic following analyses demonstrate the efficiency and simplicity of equation 
3.24 in examining the wavenumber changes of short waves riding on a current of 
longer period. Two different bottom topographies are considered: flat bottom and 
plane sloping bottom. Different current models arc applied. 

a. Waves on a Current in I l ater of Constant Depth 



i. Current as a progressive long wave 

The short wave is assumed to propagate over a flat bottom and to be 
parallel to the long wave (collincar) in the X-dircction. The current U is described by 
the orbital velocity of a progressive long period wave, 



The current will be modeled as a progressive or standing long period wave. 



cosh k ? (h+z) 



sin i ^ 



(3.25) 



U 



n 

K. 



sinh k ? h 



which corresponds to a free surface elevation, vp, = a^ sin \\i^ . Substituting into 
equation 3.24 gives 




0 



(3.26) 



equation 3.26 is easily integrated over time to give 




cos ' 4 > £ 



Constant 



(3.27) 
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where the constant is evaluated at U = 0, so that k'= constant = k s ,i.e. the 
wavenumber will not change if there is no current. Then, the general solution is 
obtained such that 



k’ 





(3.28) 



which describes the wavenumber changes due to the non-linear interaction between the 
short wave and a progressive long wave using a kinematic approach. Evaluating 
equation 3.25 at the free surface, z = 0, equation 3.28 may be written 

k' = k s ( 1 + a f k^coth k^h sin vfy ) 

It is interesting that the above result, obtained from the first principles of kincmatical 
conservation, is the same equation as (3.4 ) obtained by Longuct-I liggins and Stewart 
(1960) using a lengthy second order perturbation analysis scheme. 

ii. Current as a long standing wave 

For the Hat bottom case, equation 3.24 reduces to 



3k 1 
3 1 



+ 




(3.29) 



The current U is described by the orbital velocity of a linear standing w r ave 



U 



cosh k^ (h+z ) 
a S. U) 2.^S. sinh k^h 



sin k^x sin 



(3.30) 



which varies only in the x-dircction. Applying equations 3.29 for both x and y 
directions yields 



3k, 1 • cosh k. (h + z ) 

-TF + k s a i.Vt sinh k jh — cos k i* sin “A ’ 0 (3-3D 
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3 t 



Integrating over time gives 



cosh k (h+z ) 

~ sinh k.h Y Constant 

(3.32) 

I k ' = k = Constant 

y y 




liquation 3.32 indicates that the wavenumber component alongshore is always constant 

since no longshore component of the current exists. 

The constant in (3.32) can be evaluated at any node where = 0. At 

these locations, it is assumed that the wavenumber docs not change, i.c. k' = k. . 

x b x 

Substituting into equation 3.32, a general solution is obtained: 



I k’ = k {1 + a 0 k„coth k„h cos k n x cos w.t} , 
x s Z It Z Z Z 

x (3.33) 

) k 1 = k = Constant . 

y sy 

Therefore the wavenumber component in the x direction is always modulated in time 
and space, while the longshore wavenumber component is always constant. As a 
result, the wavenumber vector changes in time and space; i.c. the waves arc refracted. 

It is clear that if the waves arc parallel to the current, i.c. collinear in 

x-dircction, then ky = 0 and the solution reduces to the first equation in 3.33 . 

b. 1 1 a rev on Current in IJ liter of Variable Depth 

A reasonable model to describe the current on a sloping bottom is the 
orbital velocity of either a cross-shore long standing wave or edge waves. Collinear 
and oblique incident waves on a current arc included. Let us assume a plane sloping 
bottom for simplicity of analysis. 

i. Current as a cross-shore long standing wave 

On a sloping bottom, the horizontal velocity of a standing wave is 

described by equation 2.17 . Recalling the equation for wave density conservation 

(3.24), 
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(3.34) 
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8k^ 



ax. 



which requires specifying d U/d x and d 
becomes 



k„ (8 X: . The gradient of the current velocity 
n 1 



3U 

8x 



-{ 



U 



3r H 

at 



(3.35) 



Using the dispersion relationship in shallow water to differentiate the wavenumber 
components gives 





tan a 



da 

dx 




(3.36) 



and 8 ky / d y = 0, since the dispersion relationship is independent of y as the depth 
contours arc assumed straight and parallel. Applying Snell's law to describe the 
change in the angle of approach with x, d a / 8 x gives 



da 

dx 



2 ft tan a tan B 



(3.37) 



Substituting into 3.36, gives 



8k 



sx 



8x 




( sec 




(3.38) 



Substituting into 3.34 using 3.35 and 3.38 



3k x 
8 t 



+ k 



f _ U _ 1 3n f. 
s x x h 8 1 



k sec a 

1 - — i 



0 



(3.39) 
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from which 8 ky \8 y = 0 and, therefore, k'y = ky = constant. Integrating equation 
3.39 over time gives 



k' - 

x 



^ u dtI i + 5®§?0] - 



sx 



= Constant 



(3.40) 



since 

J U dt = -(l/0)~ ) 8 {]/ 8 t + constant 
Therefore, equation 3.40 becomes 



k’ 



x 



k sx 3U ri sec 2 x 

2 3t 1 2 
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sx 




Constant 



( 3 . 41 ) 



The constant in 3.41 can be evaluated at x = where current vanishes, i.c. 1 ) ( j = 0 
and IJ = 0, which gives k' x = constant = k $ . Then, for a steady oblique wave 
train, the general solution is given by 



k’ 

x 



= k 



[l + -jr 
s n 

x 



Xu, 



(1 + 



sec_a ) ouj 



9U, 

at' 



(3.42) 



If the waves and the current arc collinear (a = 0), equation 3.42 reduces to 



( 



k ' 



x 
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2 



i au, 

2 at J 

XU. 



( k 1 = k = Constant 
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This equation shows the dependence of the wavenumber changes on the surface 
elevation and on the acceleration of the current. 
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ii. Current as an Hdge Wave 

In the 3-1) model the long wave is described by an edge wave. A zero 
mode edge wave is used to simplify the analysis where the Lagucrrc polynomial equals 
unity. The horizontal velocity components are given as 



U 



a t gk yn 



-k X 

e Y sin (k y 
y £ 






V 






w , 




cos ( k y 

v 



u tt> 



(3.44) 



(3.45) 



A steady short wave train riding on an edge wave will satisfy the wave density 
conservation equation (3.24 ). Substituting into equation 3.24, using equations 3.38 
and 3.39, gives 
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Integrating both 3.46 and 3.47 with respect to time gives 
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The constant in 3.48, can be evaluated when x oo where both (J and V go to zero 
and k' x = k s . This leads to 
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Since edge waves arc progressive in the alongshore direction, equation 3.50 is similar to 
3.28 where the current was modeled as a progressive wave. The changes in the 
wavenumber and direction of short waves riding on zero mode edge waves are 
demonstrated in Fig 3.4 . The modulation in the wavenumber increases onshore 
where the wavelength goes to zero. Increasing the angle of wave incidence is shown 
to increase the wavenumber modulation. 

C. CHANGES IN WAVE ENERGY AND MOMENTUM FLUX 
1. Background 

A wave train composed of a group of individual waves propagates with the 
group velocity at which the total energy is transmitted. Energy is usually expressed in 
terms of average energy over a complete wavelength per unit surface area. The total 
energy of any wave system is the sum of both the potential and kinetic energy. 

The average potential energy lip per unit surface area is given by 
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The average kinetic energy E^ per unit surface area is estimated by 
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Tor linear progressive surface gravity waves, the average total energy E-i- is 
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Figure 3.4 Spatial changes in the wavenumber of short period waves 
due to the presence of zero mode edge wave. 
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(3.53) 
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which is equally partitioned between the potential and the kinetic energy. 

a. First Order Sea-Swell Energy Modulation 

Short period waves riding on much longer wave can be examined using the 
principle of superposition. First order modulations arc examined using the linear 
description of both waves. The linear theory is recognized as a reasonable 
approximation for its successful use in describing many observed phenomena even 
within the surf zone, (Bowen ct al., 1968; Guza and Thornton, 1980). The total surface 
elevation may then be expressed as 



= n s + 



(3.5J) 



where the subscripts s and 1 refer to sea-swell waves and long period waves 
(infragravity waves) respectively. Inside the surf zone, short waves break, dissipating 
their energy onshore. The depth-controlled wave breaking model is adopted (equation 
2.11) to describe first order short wave dissipation with zero amplitude at the 
intersection of the mean water level (MWL) with the beach. The depth in equation 
2.11 is modified to include the slowly varying depth changes due to the presence of the 
long wave, 



a = -y (h + n„ ) (3.55) 

where V|£ is the surface elevation of low frequency wave motion on a sloping beach and 
is described by either a long standing wave or an edge wave. The objective of this 
approach is to examine a first order energy perturbation (second order in amplitude) 
due to the presence of long waves. In the 2-D model, the long wave is modeled as a 
long standing wave. The potential energy for the combined short and long waves is 
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calculated by substituting the combined surface elevation given by equation 3.5-4 into 
equation 3.51 . The kinetic energy (equation 3.52) is similarly calculated using the 
total horizontal velocity, 
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where 
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(3.57) 



and Uj is given by 2.17 . It is assumed that the vertical velocities in shallow water arc 
small compared with the horizontal velocities, i.e., w < < u . The total energy Ii-p 
averaged over short and long wave periods is the sum of both the potential and kinetic 
energy and is given by 



The total averaged energy in equation 3.58 is represented by four terms at frequency 
bands co s ,o)j , (o> s + on ) and (o> s - C0j ). The lirst term in right hand side is identified 
as the short wave contribution if the short waves were considered alone. Similarly, the 
second term represents the contribution from the long waves alone. The third term is 
due to the nonlinear interaction between the short wave and the long wave and is 
represented by two side bands at the sum and difference between short and long wave 
frequencies. 
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b. Model Parameters and Preliminary Results 

Energy calculations for the combined short and long waves requires 
specifying the two parameters y and aj in equation 3.58 . The breaker index y , which 
relates the breaking wave amplitude to the local water depth, is taken to be .20 when 
using the root mean square wave height (Thornton and Guza, 1982). The second 
parameter, the long standing wave amplitude, is defined as the vertical excursion of the 
low frequency swash oscillation at the shoreline. Earlier work by Munk (1949) and 
Tucker (1950) showed the infragravity wave height to be 10 percent of the incident 
wave height, with both quantities measured in 15 m water depth. Goda (1975) 
measured infragravity wave heights in 1 m water depth and found them to be 20-40 
percent of the olfshorc incident wave height. Guza and Thornton (1982) measured 
infragravity run-up height in Torrey Pines to be 70 percent of the wind wave heights 
measured in 10 m water depth. Holman and Bowen (19S4) measured the significant 
runup height, defined as the significant vertical excursion of water level at the 
shoreline, and found it to be 60 percent of the incident significant wave height at the 
breakerline. Guza and Thornton (1985a) analyzed data acquired from three different 
experiments and found that the significant vertical swash excursion on the average was 
equal to the incident significant wind wave height measured in 10 m depth 
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(3.59) 



Therefore, the amplitude of the long wave a^is taken as .50 11^ . This result will be 
used here since equation 3.59 is weighted over three different beaches and incorporates 
the Santa Barbara data, which arc considered in the present analysis. 

Data acquired at the Leadbetter Beach experiment in Santa Barbara during 
1980 arc utilized to specify the parameters of the energy model. The best (it to the 
beach profile assuming a plane sloping beach is found to be 1:25 with a surf zone width 
of 50 m. The characteristic periods for both the short and long waves as determined 
from the spectral analysis of current meter records arc 14 and 85 seconds respectively. 
These arc typical values for the data acquired on February 4 (narrow banded waves). 
The model results arc shown in Fig 3.5 . The potential energy for long waves is shown 
to reach a maximum at the shoreline with a magnitude comparable to the maximum 
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short wave potential energy at the brcakcrline (Fig 3.5 a). The kinetic energy for both 
short and long waves decays shoreward (Fig 3.5 b). The total averaged energy decays 
onshore, with short waves dissipating all their energy at the shoreline where the long 
waves have their maximum energy (Fig3.5 c). 'flic two side band components show a 
steady growth (order of 15 percent) in total energy in the vicinity of the shoreline. 

2. Mass and Momentum Fluxes 
a. Background 

The changes in the mass and the momentum fluxes arc required to specify 
the driving forces for predicting the surf zone dynamics. Conservation equations 
developed by Phillips (1966) arc followed, in which the terms representing the mean 
and fluctuating quantities have been separated. These equations arc applicable to 
wave motion as well as general turbulent motion. No restrictions arc placed on the 
wave slopes and/or amplitudes through the use of these equations. 

The conservation of total mass per unit area can be expressed 
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(3.60) 



where i,j refer to horizontal coordinates, D is the total averaged depth of water, which 
may include the wave setup, and Mj is the total mass flux defined as 

M. / p u i dz , i = 1/2 (3.61) 

1 -h 

where Uj is the total horizontal velocity that contains a mean quantity and a 
fluctuating quantity due to the waves. The overbar denotes , at this stage, averaging in 
time over the short wave periods, but not so long as to exclude the long waves. 

fhe momentum flux equations arc derived from the horizontal momentum 
equations integrated over depth, averaged over time and using the free surface and 
bottom boundary conditions: 
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Mgurc 3.5 Theoretical wave model to predict energy inside the surfzonc, 
(a)potcntiai energy (b)kinctic energy (c)total energy. 
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(3.62) 
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The first term on the right hand side is the cltangc in the mean hydrostatic pressure 
due to the slope of the mean free water surface, and R- is the mean shear stress term. 
The first term on the left hand side is the local change in horizontal momentum flux; 
the second term contains the mean momentum flux U- Mj and the excess momentum 
due to the wave presence S" , which was termed the "radiation stress" tensor by 
Longuct-1 liggins and Stewart (1961, 1962). The advantage of using the "radiation 
stress" technique to solve physical problems is that the second order effects are 
obtained using first order wave theory, and therefore some non-linear wave properties, 
such as setup and surf beat, can be explained easily. The radiation stress for a 
progressive linear waves is given by 

h- — 1 2 M i M i 

S • • = f (pu.u. + P 6 . ■ ) dz “ Trp gD S . . — ^ 

13 _ h ; 13 13 ^ 13 p D (3.63) 



where 5 is the Kroncckcr delta. The third term can be neglected since it is of higher 
order in deep water and has a negligible contribution in shallow water. Evaluating all 
terms in equation 3.63, the excess momentum tensor for a linear progressive wave is 
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(3.64) 



where k x ,ky arc the wavenumber components and arc used in the expression to 
indicate the directions k x /k = cos a and ky /k = sin a . 

These expressions are also a good first approximation for slowly varying 
depth if a and k are interpreted as the local values corresponding to the local depth h 
(Mci, 1983). In shallow water, the group velocity c equals the phase speed c and the 

O 

radiation stress components are approximated by 
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(3.65) 
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As expected, for a standing wave on a flat bottom, the radiation stress components are 
exactly twice the values of a progressive wave (Longuct-lliggins and Stewart, 1960). 

b. Changes in The Radiation Stresses Due to the Long U'aves 

The objective now is to rederive the excess momentum due to combined 
short and long waves to develop explicit formulae. The short waves arc assumed to 
arrive obliquely at an angle u , measured counter-clockwise between the wave ray and 
the positive X-axis. The long waves arc modeled as a normally incident standing wave 
on a plane sloping bottom. The procedure of calculating the momentum flux 
components due to the presence of both the short and long waves is to substitute into 
the general formulation 3.63, which is applicable to all kinds of steady and unsteady 
flow. In the following development, the terms second order in amplitude (first order 
in energy and momentum) are retained, and all the higher order terms are neglected. 
Then, the resulting expressions arc averaged over the long wave period to eliminate the 
time dependanee of the obtained formulae. 

The radiation stress component S xx .defined as the flux of onshore directed 
momentum, is given from equation 3.63 as 
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I'hc first term under the integral is the product of the total horizontal wave velocities, 
which can be obtained by summing equations 2.6 and 2.17 . The limits of integration 
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arc taken from the bottom to the still water level (SWL) at l| = 0, which is consistent 
with the analysis being second order in wave slope. Integrating over depth and 
averaging in time over long wave periods yields 
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(3.69) 



The pressure term in equation 3.68 is obtained by integrating the vertical momentum 
equation (neglecting the vertical stresses) from any depth z to the free surface (see 
Phillips, 1966), which gives 

n g n g n 2 

P| = p{ / gdz + ^ / wdz + / uwdz - pw T } (3 70) 

z z z j z 



The first term is the hydrostatic pressure. Integrating from the bottom to the free 
surface and averaging over short and long wave periods gives 
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The second term is similarly evaluated, but it was found to vanish when averaged over 
the long wave period. The third and fourth terms give 
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The time averaged radiation stress S xx is obtained by adding equations 
3.69, 3.71, 3.72 and 3.73 : 
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(3.74) 



The first term on the right hand side is identified as the short wave contribution if 
considered alone, which is in agreement with equation 3.64 obtained by 
Longuct-I Iiggins and Stewart (1962). The second and the third terms arc the 
contributions by the long waves and the side bands. 

The transverse component of radiation stress, S Y is defined as the onshore 
flux of longshore momentum across a vertical plane any distance offshore. This 
component can be calculated from the direct evaluation of the momentum flux as the 
product of the total horizontal velocities, u-p and Vy , 



S y X " _J U T V T dZ (3.75) 



A cross-shore long wave is considered for simplicity, i.c. no longshore velocity 
component, while an oblique short wave is considered. Then, equation 3.75 becomes 
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The transverse component Sy X is the driving force for the longshore current inside the 
surf zone. When considering the cross-shore long waves, the primary change in Sy X 
arises from the amplitude modulation of the short waves, which is a function of the 
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long waves. Integrating equation 3.76 over depth and averaging over the short wave 
period gives the time dependent form, 
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Averaging over the long wave period to eliminate the time dependency, the transverse 
component Sy X inside the surf zone is 

c 1 r 2, 2 ^ 1 2 2 t 2 . . , . (3.78) 

S yx = 2 pg Y h + 2 Y Vo (x)}sin a cos a 

The first term on the right hand side is in agreement with equation 3.67, while the 
second term can be identified as side bands. In equation 3.78 , it is noted that the 
cross-shore long wave by itself does not contribute to S , but their interaction with 
the short wave results in the side bands that increase the Sy X term. 
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IV. DATA ANALYSIS 



The analysis of field data is performed to test the hypotheses that the nearshore 
dynamics are due to both the short and long waves. Outside the surf zone, waves 
arriving in groups cause a change in MWL (Fig. 4.1 ) and drive a forced long wave 
(Longuet-Higgins and Stewart, 1962). At the breaker zone, the short waves break and 
their heights become depth dependent, obeying the saturation curve. The water depth 
is modified by the presence of the long waves. Therefore, the short waves are expected 
to be modulated by the infragravity waves inside the surf zone. The experiments are 
briefly discussed, followed by analyzed results. 



A. DESCRIPTION OF THE FIELD SITES 

Two experiments were performed as a part of the Nearshore Sediment Transport 
Study (NSTS). The objective of NSTS was to develop an improved engineering 
formula to predict sediment transport on beaches. Field sites were selected for their 
simple topography of straight and parallel depth contours, which match the 
assumptions often used for nearshore dynamics. The first experiment was located at 
Torrey Pines Beach (T.P.), San Diego, California and conducted for one month in 
November, 197S. Torrey Pines is a gently sloping beach (tan p = 0.002) with nearly 
straight and parallel contours. The bottom is composed of moderately sorted, 
fine-grain sand (mean diameter 0.15-0.2 mm). The waves break as spilling or mixed 
spilling-plunging. A wide variety of wave and weather conditions were encountered 
during the experiment, from small to large (2 m) waves, from very narrow-band swell 
to wide-band sea, from calm to windy days (> 10 m/s). A total of 42 sensors were 
deployed over an area 520 m long parallel to the shore and up to 500 m offshore. 
Surface elevation and horizontal, orthogonal velocity components were measured along 
a shore-normal transect from offshore at the 10 m depth contour to the shore. 

Leadbettcr Beach, Santa Barbara (S.B.), California was the second NSTS 
experimental site (February', 19S0). Leadbettcr Beach has relatively straight and 
parallel nearshore depth contours with well-sorted fine to medium size sand. The mean 
nearshore slope varied between 0.017 and 0.05 during the experiment, depending on the 
wave climate. Well-developed cusps occurred during the beginning of the experiment. 
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Figure 4.1 Schematic diauram of. short and long wave models outside 

and" inside the surf zone. 
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No major offshore bar system was apparent. The shoreline has an unusual east-west 
orientation along a predominantly north-south coast. The open ocean waves are 
limited to a narrow window of approach (± 9° centered on 249°) because of 
protection from Point Conception to the north and the Channel Islands to the south. 
The generally highly filtered ocean swell type waves from almost due west must make a 
right angle turn to approach the beach normally. The result is often a narrow band (in 
frequency and direction) spectrum of waves approaching at large oblique angles to the 
bottom contours in the surf zone and driving a strong longshore current. Rod and 
level surveys of the beach profile were obtained daily. An array of current meters and 
pressure sensors were deployed extensively in a similar manner as at T.P. Severe 
storms occurred during the experiment, causing significant changes in the shoreline 
configuration and beach profile. Despite the destruction of most of the instruments 
during the storm, a very substantial and valuable data was recorded (Gable, 1980). 

Data were acquired at a sample rate of 64 Hz for several hours each day, then 
low-pass filtered and reduced to 2 Hz. The sensors, data acquisition system and 
experimental sites are fully described by Gable (1979, 1980). 



B. DATA ANALYSIS 

Data acquired from T.P. and S.B. have been extensively analyzed by different 
investigators. Incident W’ave analyses include shoaling, setup, swash oscillation, 
runup, velocity moments, longshore current and sediment transport. Infragravity 
waves were studied for runup, identifying both the cross-shore and longshore structures 
and type of the dominant motion, whether they are leaky and/or edge waves (Guza and 
Thornton, 1982, 1985a, b; Oltman-Shay and Guza, 1986). As mentioned before, the 
dynamics of the infragravity waves were not previously studied. Furthermore, their 
effect on the incident waves, especially inside the surf zone, is of interest to understand 
the nearshore processes. It is hypothesized in this study that the interaction between 
the short and infragravity waves generates other waves (side bands) at the sum and 
difference frequencies. These waves occur at the ncigbourhood of the short wave 
frequency, and they arc of importance in driving longshore current and sediment 
transport at the shoreline. The existence of both the infragravity waves and the side 
bands arc investigated utilizing records of wave staffs, current meters and pressure 
sensors. Coupling between the incident waves and the infragravity waves arc tested via 
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spectral and cross-correlation techniques, and different theories arc discussed. Since 
the present developed model is monochromatic, only days of narrow banded waves arc 
considered. Three days are selected from each experiment (total off) days): 3, 4, and (> 
February 19S0 at S.B. and 10, 20, and 21 November 1978 at T.P. 

1. Spectral Analysis 

Velocity and pressure time scries of 68-minutc length arc used to compute 
spectra with 32 degrees of freedom and a resolution of 0.004 Iiz. The power spectra 
arc computed for all the available working sensors at each selected day. The spectra 
arc plotted on linear scale to demonstrate the relative importance of each peak and 
normalized by the variance (area under spectra). In Fig 4.2 , a well defined swell peak 
is seen at a frequency around f $ = .06 Iiz in most of the instrument locations during 
10 Nov., T.P., and in Figures (4.3, 4.4) at frequency f s ~ .07 Iiz in all the analyzed 
instruments in 3 and 4 Feb., S.B. A consistent valley near 0.05 Iiz is present in almost 
all spectra, separating the incident wave band from the low frequency band. Within 
the infragravity band, a pronounced low frequency peak is shown in most of the 
spectra at fg = .017 Iiz for 4 Feb, S.B. and .0195 Iiz in 10 Nov., T.P. 

'I he infragravity energy has been studied extensively since the first 
observations by Munk (1949) and Tucker (1950). Huntley ct al. (1981) identified a 
progressive low mode edge wave in the longshore direction that satisfies the dispersion 
relationship of edge waves. Oltman-Shay and Guza (1986) extended the work of 
Huntley et al. (1981) and were able to determine the energy contained in each edge 
wave mode. The low mode edge waves (n ^ 2) were found to dominate the longshore 
current energy while leaky mode or high mode edge waves dominated the cross-shore 
direction. They attributed the daily changes in infragravity energy levels to the 
variation in incident wave conditions and the resonant forcing of edge waves described 
by Gallagher (1971). Guza and Thornton (1985a) observed surf beat in 3 field sites 
including S.B. and T.P. They found a phase shift of n / 2 between the surface elevation 
and the horizontal velocity within low frequency motions (f ^ .05 Iiz) which indicates 
the predominance of the standing wave structure of surf beat. This result agrees with 
Suluiyda (1974). Guza and Thornton (1985a) rejected the hypothesis that surf beat is 
primarily generated by bore-bore capture within the surf zone (as hypothesized by 
Bradshaw, 1980 and others) because bore-bore capture models do not predict incoming 
surf beat energy outside the surf zone and do not describe the gradual increase in the 
infragravity band onshore. 
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Figure 4.2. Normalized power spectra at different locations in 10 N'ov., T .1*., 
showing an onshore relative increase in the infragravity energy level. 
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The transformation of the incident wave spectra across the surf zone is shown 
in Figures 4.2 and 4.3 . A dramatic shoreward increase in the low frequency band 
energy is accompanied with rapid attenuation in the progressive incident wave energy 
due to breaking as the depth decreases. The infragravity energy can exceed the wind 
wave energy by a factor of 2.5 in the inner surf zone (Wright ct al. , 1982). In storm 
conditions, Holman ct al. (197S) observed a considerable increase in the infragravity 
energy levels responding to the changes in incident wave conditions. It is clear from 
Figures (4.2, 4.3 and 4.4) that the spectra closer to the shoreline is dominated by the 
low frequency band. 

The theory of wave-wave interaction presented in Chapter 3, by which side 
band waves arc generated as a result of the non-linear interaction, can be tested using 
the available data, lilgar and Guza (1985), on studying the nonlinear dynamics of the 
shoaling waves using bispcctral analysis, showed that the low frequency motion at the 
beach face is significantly non-lincarly coupled to higher frequency modes seaward the 
surf zone around the peak of the power spectra. They did not resolve possible side 
band generation in their analysis, choosing instead to obtain higher stability (low 
resolution). Using higher resolution here, peaks at frequencies (f $ ± fj, ) are found in 
Fig. 4.5 which are identified as upper and lower side bands. This supports the 
hypothesis of the short wave modulation derived in equation 3.11 . Evidence of side 
bands is present for almost all the measurements on all days analyzed, but usually only 
the upper side band at frequency ((' s + Q ) is present (see Figures 4.2, 4.3 and 4.4). 
The side bands arc not always statistically significant, however, since at least the upper 
side band is present for all narrow band days, lienee, it is concluded that the side 
bands arc real. It is not clear why the lower side band is absent. The suppression of 
the side band may be strong lower frequency breaking or it may be caused by the 
seaward radiation of energy by the outgoing component of the standing long waves. 
FIgar (1985) showed that the low frequency reflected waves will have phase 
relationships relative to their incoming parts that tend to decrease bispcctral level. 
Guza ct al. (1984) show an increase in coherence between the incident wave envelope 
and the long standing waves when removing the reflected wave. 

The time history of energy transferred between swell peak and the upper side 
band is demonstrated in Fig. 4.6 using almost a four-hour time scries, 'flic long record 
is divided into 64 segments each of 17.1 minutes with a 50% overlap to calculate the 
spectra. Fnergy appears to be transferred from the incident wave peak to higher 
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Figure 4.3 Normalized power spectra at different locations 

in 3 I cb., S M. 



67 



NORMALIZED ENERGY DENSITY 

0.0 20.0 10.0 C.O 20.0 10.0 0.0 20.0 




Figure 4.4 Normalized power spectra at dilTerent locations 

in 4 Feb./S.U. 
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Figure 4. .5, Normalized spectra showing significant peaks at 
frequencies of dominant long wave fjj, short wave l‘ s and side bands. 
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frequency components. This leads to a faster growth of the upper side band than for 
the lower side band. The incident wave peak drops to a minimum, strengthening the 
upper side band. A similar mechanism could happen if lower side bands exist, which 
may transfer energy to low frequency components. 

Using the same data set, Guza and Thornton (l9S5b) computed the velocity 
variance in 20 Nov., T. I’, at each instrument location across the surf zone Fig. 4.7 . 
Their result show no substantial differences between the observed and the predicted 
values when summing the potential and kinetic energy for short waves riding on long 
standing waves (Fig. 3.5). The decrease in short wave energy is almost balanced by the 
increase in low frequency energy resulting in an almost constant cross-shore velocity 
variance (Guza and Thornton, 1985b). Wright ct al. (1982) obtained a similar variance 
distribution on Australian beach. Although the spectral shape of 3 Feb is similar to 4 
Feb, waves in the 4 Feb arc more energetic (1L:„ = 80 cm). Therefore, strong 
non-linearities arc expected. 

2. Cross-Correlation Analysis 

The objective of this analysis is to cross-correlate the incident wave envelope 
with the long waves to test the hypotheses that long waves are forced outside the surf 
zone and short waves modulated inside the surf zone. The cross-correlation function, 
Cj^ , is a measure of the degree of linear relationship (association) between two data 
sets in the time domain. Tucker (1950) correlated the incoming wave envelope and the 
corresponding long wave observed by a bottom mounted pressure sensor at a distance 
800 nr offshore. He obtained a negative maximum correlation corresponding to a 5 
minute time lag. The time lag associated with the maximum correlation was 
approximately the time required for the incident wave group to reach the shore with 
velocity Cg and for the associated long wave to travel back with phase velocity J gh. 
The negative correlation was described by Longuct-Higgins and Stewart (1962) as a 
low frequency forced wave response due to the incident wave groupness that has a 
phase shift of it with the envelope of the high frequency waves. Hence, under high 
waves, a long wave trough is expected and vice-versa. 

lluntly and Kim (19S4) obtained a positive correlation between long waves 
and the short wave envelope at zero time lag from current meter measurements taken 
outside the breaker line. Their results agree with the forced wave model of 
Longuet-1 Iiggins and Stewart (1962). For velocity measurements, a positive correlation 
is expected with the x-dircction positive offshore. 
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Figure 4.6 Spectral time series of current meter OODx with 50% overlapping 
calculated from four hour record, (4 Feb, 1980, Santa Barbara). 
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Guza, ct al. (1984) obtained a negative correlation outside the surf zone and a 
positive correlation inside the surf zone when comparing the surface elevation of the 
short wave envelope with long waves. No explanation was given for change in signs of 
the correlation between outside and inside the surf zone. 

Surface elevation time series are inferred from current meter data from 
experiments at T.I\ and S.B. using the complex Fourier spectra of the horizontal 
velocity component, U(f). Current meter records were first linearly detrended to 
exclude the effects of the rising and falling of the tides, then high-pass filtered using a 
cut-oil' frequency of 0.05 Hz to exclude surf beat, followed by a low-pass filter with a 
high frequency cutoff (0.31 Iz) to exclude capillary waves. The filtering was 
accomplished by Fourier transforming the signals and zeroing the Fourier amplitude 
coefficients in the filtered frequencies. The complex surface elevation spectrum, X(f), 
is calculated applying the linear wave theory transfer function, llj-(f) 

X ( f ) = H f (f ) U (f ) 

f (4.1) 

Then, the complex surface elevation spectrum is inverse transformed to obtain the 
surface elevation time series. The entire 68-minute-rccord was transformed at one time 
to minimize the end effects, which result in spectral leakage, and to obtain maximum 
resolution for very sharp roll-off at the filter cut-offs. Surface elevations were also 
computed from pressure signals in a similar manner by transforming the pressure 
records using linear theory. 

Using the surface elevation time series obtained in this manner, the correlation 
between the incident wave envelope and the corresponding low frequency motion is 
calculated. The low frequency motion is obtained by low-pass filtering the surface 
elevation time scries. The high frequency wave signal is the residual after subtracting 
off the low frequency energy. The short wave envelope is obtained by demeaning and 
squaring the high frequency signal, followed by the same low-pass filter used before. 
Segments of the time series for the wave envelope and the associated low frequency 
motion suggest positive correlation for the wave sensor inside the surf zone (Fig. 4.8) 
and negative correlation for the current meter derived signal outside the surf zone (Fig. 
4.9). 

The cross correlation function is calculated using the following formula 
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Figure 4.7 Spatial distribution of band passed velocity variance in 20 Nov., 1978. 
b 1 Low is (0 < 1 < 0.05 I Iz.) and high is (0.05 <f< .5 1 Iz) 

(alter Guza and Ihornton, 19b5b). 
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Figure 4.8 Segment of the time series for the wave envelope (solid curve) and the 
associated low frequency (dotted curve) for wave stalF W38 (inside the surf zone) 

in 21 Nov., 1978, Torrey Pines. 
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Figure 4.9 Segment of the time scries for the wave envelope (solid curve) and the 
associated low frequency (dotted curve') for current meter CO/x (outside the surf zone) 

in 4 Feb., 1980, Santa Barbara. 
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where t = k t is the time lag, t is the sampling time interval, n is 
number of data points, <T| <r e arc the standard deviations for the low frequency signal 
and the envelope, and x, y are the means for each of the time series. The time lag is 
positive when the long waves arc in advance of the wave envelope, and it is negative 
when long waves lag the wave envelope. 

A test of significance (null hypothesis) based on the "Fishcr-Z 
Transformation" is performed with the statistic z = J n-3 /2 In 0 + C R / »- Cr ) 
(Miller and Freund, 1977). At level of significance a = 0.05, the calculated z always 
exceeds z tf /2 (= 1.96), since the used record is long (S192 data points). Thus, the null 
hypothesis is rejected and the correlation is significantly different than zero. 

The correlation function for the wave staff \V3S (T.P.) is show n in Fig. 4.10 . 
The maximum correlation is +0.49, corresponding to a time lag of -4.0 sec. The 
function shows a strong periodicity at about 40 see. The correlation coefficient is 
calculated with 95% confidence to be significantly different from zero. Fig. 4.11 shows 
a negative correlation for C07x (S.IL), indicating the characteristics of a forced wave 
response outside the surf zone. 

The maximum correlation coefficients and the corresponding time lags at each 
instrument location in 10 and 21 November, T.P. arc given in Table 1 and 11 . 
Similarly, tables 111 and IV arc given for 4 and 6 Fcb.,S.B. Lower values of imply 
that the two random variables may be independent, that there is no linear relationship 
between the two signals, suggesting a free long wave (Iluntly and Kim, 1984) or the 
possibility of some non-linear relationship. The low values of the negative Cn can be 
attributed to the fact that the n phase shift is only a deep water ease, while the deepest 
instrument used is 3 m. The other reason could be the significance of the phase-locked 
reflected low frequency waves decreasing the coherence even for low values of 
reflection coefficient (lilgar, 19S5). Hlgar and Guza (1985) indicated that the 
infragravity modes do not appear to be bound with a fixed phase relationship to the 
high frequency wave groups, since thci r biphase evolves towards lower values as waves 
shoal. 
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I'igurc 4.10 The cross-correlation function for W38, 21 Nov., 1978, 
Torrey Pines. A positive maximum correlation (.49) corresponds to -4.0 time lag. 
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Figure 4.1 1 The cross : correlation function for C07x, 4 Feb., 1980, 

Santa Barbara. A negative maximum correlation (.29) corresponds to +3.0 time lag. 
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It is of interest to observe a consistent pattern of C|^ , being negative outside 
the surf zone and positive closer to the shoreline (Figures 4.12 and 4.13). '1 he 

statistically significant values of in the nearshore zone at zero time lag indicated 

significant coupling between the high and low frequency motions. A zero correlation 
in big. 4.12 is found to coincide with the approximate location of the breaker line. The 
zero correlation could be due to the two signals being non-linearity related, 
randomness (noise; introduced into the wave signal by the breaking wave, a possibility 
that the breakpoint could be at a nodal point for the long standing wave, independence 
of the two signals, or a transition location between forced and free long waves. Fig. 
4.13 shows a consistent negative correlation outside the surf zone and a positive 
correlation shoreward of the midsurf zone. This could mean that the wave breaking 
docs not destroy completely the wave groupness until some location inside the surf 
zone. This location might be a transition between the forced and free long waves. 
Flgar (1985) showed that surf beat biphases decreases as waves shoal. On some days, 
the correlation in Fig. 4.13 shows a rapid decrease when waves get closer to the 

shoreline. This might be attributed to the bore-bore capture mechanism. Masc and 

Iwagaki (1984; illustrated that a wave (bore) could be captured by a successive bore 
before the bore reaches a maximum runup height. Therefore, the number of runup 
waves and frequency of runup will be reduced compared to that of incident waves. The 
ratio of number of waves is given using the surf similarity parameter and shows a 
decrease in number of waves reaching a gentle slope beach. 

The positive correlations obtained inside the surf zone are demonstrated by 
the simple model developed in chapter 3 by linear superpose the two waves. The 
surface elevation time scries inside the surf zone is simulated using 

n T = v(h + n ^ ) cos (4.3) 

The wave envelope and low frequency motion arc obtained as described above. A 
positive correlation is always shown at zero time lag (Fig. 4.14 ), which agrees with the 
observations. 

In the following chapter, the theoretical results obtained in Chapter 3 arc 
applied to derive wave setup, longshore currents (steady and unsteady) and sediment 
transport due to the combined short and long period waves. 
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Figure 4.12 Spatial distribution of cross-correlation coefTicients 
snowing a dominant negative correlation outside the surf /.one . 
and a positive correlation inside the si.il /one. B. h. is the I reakerline. 
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Figure 4.13 Spatial distribution of cross-correlation coefficients 
showing a dominant negative correlation outside the surf zone 
and a positive correlation alter midsurf zone. 
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(b) 




Figure 4.14 Theoretical model to examine the positive correlation 
between the wave envelope and the associated low frequency motion within the 
surf zone, (a) segment of time series (b) cross-correlation function. 
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V. APPLICATIONS 



A. WAVE SETUP AND SETDOWN 
1. Background 

In the nearshore zone, it is observed that waves during shoaling and breaking 
produce a variation in the mean water level (MWL). This variation in MWL is 
considered a primary cause of nearshore circulation, such as rip currents. The change 
in MWL, , is characterized by two regions (a) a gradual depression of the sea water 
level starting from offshore, and reaching a maximum at the brcakerlinc, and (b) an 
increase in the MWL shoreward of the breaker line. The depression is termed 
"setdown" and the upward slope is called "setup." 

To study the wave setdown and setup, the X-componcnt of the horizontal 
momentum (lux equation is considered where the waves arc assumed to be steady state 
and propagating normal to the shoreline over straight and parallel depth contours, 



Equation 5.1 indicates that the change of the excess momentum llux due to wave 
motion in the cross-shore direction is mainly balanced by the mean water slope dlj/dx, 
assuming no gradients in the Y-dircction and neglecting the frictional forces. 

It is common to separate the nearshore zone into two dynamical regions of 
"outside" and "inside" the surfzonc. 

(1) Outside the SurfZonc : 

If waves are assumed to propagate onshore without any energy loss, then the 
onshore component of the radiation stress is expected to increase steadily as the depth 
decreases. Therefore, the mean water level is lowered (set down) by the presence of the 
unbroken waves. Longuet-I liggins and Stewart (1962) integrated equation 5.1 by 
assuming that Tj” < < h to obtain 



as 



XX 




(5.1) 



n 



1 a 2 k 

2 sinh 2 kh 



(5.2) 
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This equation is second order in local wave amplitude. The setdown increases steadily 
until the breaking zone. Without referring to the radiation stress, Longuct-1 liggins 
(1967) also derived the above solution using Bernoulli's integral equation. This simple 
technique will be illustrated and applied below. 

(2) Inside the Surf Zone 

As waves break, the wave amplitude inside the surf zone is assumed to be 
controlled by the total local water depth including the wave setup q , simulating a 
spilling breaker type, i.c. 



a = y ( h + n ) 



(5.3) 



The cross-shore component of the radiation stress, S xx , is a function of the energy 
density which decreases inside the surf zone and reduces to 3/2 E in shallow waters 
such that 



S 



XX 



3 2 
q-PTY 



(h 



+ n ) 2 



(5.4) 



Both the gradients of S xx and the forcing in equation 5.1 have a negative slope. 
Therefore, the change in the MWL inside the surf zone is expected to be positive, 
representing a wave setup. Integrating equation 5.1 gives 



n = N(h b -h) + n b 



(5.5) 



? 1 

where N = (1+ 2/(3 y )) . The subscript b indicates the value at the brcakerlinc 

where = - y 1 I b /8 is determined from equation 5.2 . Equation 5.5 shows an 
increase in the mean water level with decreasing depth, balancing the decreasing 
radiation stress duc to wave dissipation inside the surf zone. It should be noted that 
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Bernoulli s integral outside the surf zone is not applicable inside the surf zone since the 
motion is no longer irrotational. 

Bowen ct al. (1968) demonstrated the above theoretical results by careful 
measurements in the laboratory. Experimental results agreed well with the theory. 
However, just outside the breaking point, their experimental values were less than the 
theoretical predicted since the theory docs not consider energy dissipation. Also, very 
close to the shoreline, Bowen et al. (1968) found a residual wave height (possibly due 
to a standing wave) to persist and corrected the above equation by letting the setup 
slope be tangential to the beach slope, i.c. dq/dx - dh/dx as (h + q)-^ 0 . Van 
Dorn (1976) found that the setup gradient at the shoreline is actually steeper than 
predicted by equation 5.5, which agrees with Bowen ct al. (196S). Battjes (1974) and 
Battjes and Janssen (1978) used a random wave model, which predicts smaller setup 
values. 

In addition to the wave setup laboratory investigations, a few field 
measurements have been reported. Dorrestein (1961) measured the change in mean 
water level across the surf zone during times when the significant offshore wave height 
ranged from 0.8 to 1.6 m. lie measured a maximum setup of 0.15 m. Guza and 
Thornton (1981) measured the setup at the shoreline when the significant offshore 
wave height varied from 0.6- 1.6 m. The maximum setup at the shoreline was found to 
be about 0.17 IIq , where I Iy is the deep water significant wave height. Holman and 
Sallcnger (1985) measured the setup under conditions of incident significant wave 
height varying from 0.4- 4.0 m. They found the setup to vary linearly with the surf 
similarity parameter Cq = P (IIq /Lq )'*' L . 

2. Changes in Setup and Setdown Due to the Presence of Long Waves 

Changes in the wave setup and setdown due to the presence of long waves can 
be demonstrated by considering the same wave field used before, i.c. oblique short 
waves riding on a normally incident long period standing wave. The wave setdown is 
evaluated applying Bernoulli's integral equation (Longuct-IIiggins, 1967). first, the 
vertical momentum equation is integrated over depth and averaged over time. 





(p +pw ) | - pgn = 0 

Z 



(5.6) 
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where p is the pressure and the overbar denotes mean quantities. Next, the time 
averaged Bernoulli's integral equation is evaluated at z = 0, assuming the flow is 
irrotational, 



, 1 2 
P | + jp(u 

z = 0 



-2 - 2 , 
+ v + w ) 



+ Constant = 0 



z = 0 



(5.7) 



The pressure term is eliminated by combining equations 5.6 and 5.7 to obtain 
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(5.8) 



The difference in mean sea level 8 i] is obtained by applying equation 5.8 at two 
different locations (xj ,y j ,0) and (X 2 ,>’2 .0), which avoids evaluating the constant. 
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(5.9) 



For the combined short and long standing waves, equation 5.9 is used to evaluate the 
setdown as follows. Hach velocity component is expressed as the sum of the short and 
long standing wave velocities as given before in chapter 2. Let location 2 be in deep 
water (h 00 ) where the setdown is zero. Substituting into equation 5.9 and 
averaging over the long wave period gives 



n = _ s " L /„ 2 _ 2 n ( 5 - 10 ) 

2 sinh 2kh 2g U Z w £ 

where subscripts s and 1 refer to short and long waves. The first term in equation 5.10 
is identified as the the setdown solution (equation 5.2) when considering the short 
wave only. At the' breakerline, the depth controlled breaking model is utilized and 
setdown reaches a maximum to give 
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(5.11) 
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where the subscript b refers to conditions at the breakcrline. 

Wave setup can be evaluated by spatially integrating the X-momentum flux 
equation, where the radiation stress component S xx for the combined wave fields is 
given by equation 3.74 . The constant of integration is evaluated by applying the 
obtained setup equation at the breakcrline where the change of the mean sea level is 
known from equation 5.11 . The resulting equation is quadratic in i) and has two 
roots. One root is disregarded since it unreasonably predicts setdown at the shoreline. 
The second root gives 



n = -(h + v^h 2 -q ) 



(5.12) 



where 



Q 



- [r 'b + 2 V'b T 2 wt '“b 



+ i(3y 2 (h 2 - h 2 ) 



+ a 2 <5i DtJo 



(5.13) 



+ 2a 2 (J 2 (X) - J 2 (X) ) } ] 



Using Longuet-1 Iiggins' setup model for short waves as a reference, the setup is 
relatively increased closer to the shoreline due to the presence of the long standing 
wave (Fig 5.1 ). However, the setdown is also increased at the brcakerline compared 
with the reference model. The relative increase in setup appears to be over-predicted 
and suggests the importance of including damping terms, such as frictional dissipation 
and/or percolation, in the swash zone. 
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Figure 5.1 Wave setup due to. both incident and infragravity waves 
compared with Longuct-1 liggins and Stewart model (1962). 
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B. LONGSHORE CURRENT 
1. Background 

When waves obliquely approach the surf zone, a mean current parallel to the 
shoreline is generated, which is important in transporting sediment alongshore. In the 
simplest formulations, the longshore current is assumed steady and two dimensional, 
being independent of the longshore topographic changes. The linear Y-momentum flux 
equation is simplified such that the gradient of the radiation stress component d S yx / dx 
is balanced by the bottom shear stress, 



ifvx = _ (5.14) 

ox y 



Outside the surf zone, the wave energy is assumed non-dissipative and the excess 
momentum flux in the Y-dircction is a conservative quantity. Therefore, no driving 
force is expected outside the surf zone for generating a longshore current (Bowen, 
I960). Inside the surf zone, energy is dissipated due to wave breaking, causing a 
change in the momentum flux. The onshore changes in momentum llux drive a 
longshore current primarily confined to the surf zone. 

Numerous models of longshore current on a long straight beach have been 
done since the introduction of the radiation stress concept. For obliquely incident 
waves on plane beaches, analytical solutions were derived for monochromatic waves by 
Bowen (1969) and Longuet-I liggins (1970 a,b) and for random waves by Thornton and 
Guza (1986). Numerical treatments (e.g. Thornton, 1970; Jonsson ct ah, 1974; 
Madsen et ah, 1978); Wu and Liu, 1984, and others) have employed arbitrary beach 
profile and/or including setup and/or non-linear inertial terms. In monochromatic 
wave models, lateral mixing is introduced to smooth out the discontinuity in the 
longshore distribution due to the intensive current shear at the breakerline (Bowen, 
1969; Thornton, 1970; Longuet-I liggins, 1970a). For random waves models, there is 
no need to include the lateral mixing (Collins, 1970; Battjcs, 1972) since the waves 
transform and break over a considerable horizontal distance resulting in a smooth 
change in the rate of energy dissipation and S yx . 

Periodic variation and unsteadiness in the longshore current measurements 
have been reported 'by Wood and Meadows (1975) and Meadows (1976), but neither 
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analytical nor numerical solutions have been attempted. Meadows (1976) suggested 
three distinct longshore velocity components, a steady component and two fluctuating 
components, with one at a short wave frequency and the other at low frequency. 
Guza and Thornton (1978) and Holman and Bowen (1984) observed low frequency 
components in the longshore velocity spectra. Since the longshore current is a 
time-averaged velocity, Guza and Thornton (1978) stated that an appropriate temporal 
averaging time for mean longshore current is unknown. The above studies suggest that 
the low frequency wave motion can be important inside the surf zone. 

Infragravity waves arc included in the following wave description inside the 
surf zone, allowing a more complete description of the wave field. 'The changes in S„ v 
due to obliqueness of the short wave and the total bed shear stress arc derived. Then, 
the steady and unsteady longshore current formulations arc considered on plane 
sloping beaches. 

2. Wave Refraction and Radiation Stress 

The depth contours arc assumed to be straight and parallel, and obliquely 
incident waves arc therefore refracted according to Snell's law, 



sing 

c 



sin a. 

= Constant 

c b 



(5.15) 



where c is the phase speed. The S v component for the combined short and long 

yx 

waves is given by equation 3.77 . Assuming a small angle of wave incidence (cos a ~ 
1), Snell's constant (sin u /c), which contains all the wave angle information, can be 
introduced into equation 3.77 to yield 
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(5.16) 



Another approach for deriving equation 5.16 can be shown by using equation 3.67, 
where 
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(5.17) 
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By assuming a small angle of wave incidence, 5.17 is reduced to E c (sin a /c) where E 
is given by 1/2 p g y2 (h + i|j ) fc . This will give the same solution as equation 5.16 . 
This result can be attributed to the fact that standing waves by themselves do not 
contribute directly to the longshore momentum flux since they are oriented in the 
cross-shore direction, and no S forcing is expected for a normal standing wave. The 
contribution of the infragravity waves manifests itself by the generation of momentum 
at the side bands. Averaging equation 5.16 over the long wave period gives 
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(5.18) 



That is to say, the S^ x forcing is increased by the formation of the side bands (second 
term in bracket). The longshore current forcing term d Sy X / d x for both steady and 
unsteady cases is compared in Eig. 5.2 with the Longuet-I liggins ( 1 97(>a) model where 
only short waves are considered. In general, the combined forcing by short and long 
waves increases in the vicinity of the shoreline due to the relative growth of the side 
bands. 

3. Bottom Shear Stress 

Inside the surf zone the bottom stress can no longer be neglected, especially 
when considering the change of momentum flux in the Y-direction. Assuming the 
quadratic shear stress formula, the averaged alongshore bottom shear stress, exerted by 
the oscillating wave motion is given by 



x 

y 





(5.19) 



where Cp is the friction coefficient and u-p is the vector sum of the current and the 
in viscid orbital velocities for both short and long waves measured just above the 
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Figure 5.2 Gradient of the longshore current forcing across the 

surf zone. 
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bottom boundary layer. The absolute value is necessary to insure that the stress 
reverses with velocity reversals. Longuet-l liggins (1970a) assumed a small angle of 
wave incidence and a weak current compared to the water particle motion of the 
waves, and he was able to simplify 5.19 to 
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p c f |u T | v 



(5.20) 



where V is the mean motion in the longshore direction which is assumed uniform over 
depth. The weak current assumption has been shown to be a reasonable 
approximation for many field applications (Liu and Dalrymple, 1978). The total 
velocity vector u-j- can be resolved into its components in the cross-shore and 
alongshore directions: 
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(5.21) 



where u $ and Uj are the X-component of the orbital velocities and i,j are the unit 
vectors in the x and y directions respectively. U and V are mean velocities. Assuming 
uniform currents alongshore, the mean cross-shore velocity U is zero according to the 
conservation of mass, since the beach forms a boundary in the X-direction. 
Considering an oblicjnc short wave riding on a cross-shore long standing wave, i.c. Vj 
= 0, equation 5.21 can be rewritten 



u„ 



= { ( u cos a + 
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V}J 



(5.22) 



The absolute value of the total velocity vector is 
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Assuming a small angle of wave incidence (« < 10° ), cos a can be approximated to be 
unity, which simplifies 5.23 to 
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(u 
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2u V sin 
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(5-2 A) 



It follows that u $ sin a is small and can be neglected, then 



l“ T l = V (U S + .“I 1 " + y2 

which can be rearranged into 

|u T l = |u s + u 4 | V 1 + 



(5.25) 



(5.26) 



Note that u $ + uj = Uy i according to equation 3.56 . Then, using the weak current 
assumption, V/uy < < 1, to simplify the analysis gives 




u s + u t l 



(5.27) 



Since the two velocities u $ and iij are almost collincar (both acting almost in the 
X-direction), equation 5.27 gives the basic relationship 




Equation 5.28 can be evaluated using 3.57 and 2.17 and averaged over the short wave 
period to give 
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(5.29) 
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